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Introduction 


•  High-order  in  space  is  now  commonplace 

•  High-order  in  time...  not  so  much... 

*  Is  this  sufficient?  Is  high-order  in  time  needed? 

*  Limiting  Fact  There  are  no  A-stable  backward-difference 
formula  (BDF)  methods  with  >  2nd  -order  accuracy 

*  Thus,  multistage  methods,  like  Runge-Kutta  (RK) 
methods,  must  be  used  for  3rd-  and  higher-order 

•  Explicit  RK  methods  are  not  amenable  to  stiff  problems 

Objective:  T o  find  optimal  diagonally-implicit 
Runge-Kutta  time  integrators  for  use  with 
high-order  spatial  discretizations 
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ife#  Governing  Equations 


Dual  Time  Stepping: 

<9Q  dQ  OF  i  dV 


dr+^+d^-^+H  Q  =  iP  PUi  pe  o] 
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F i  =  \pui  pUiUj+pSij  Uiph0]  where  h0  =  e0  + 


T 


Quasi-linear  Form: 
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A  =  diag  {ui  +  c,Ui,Ui  —  c } 


Residual  Form: 
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ife#  Spatial  Discretizations 


Central  Differences  with  added  artificial  dissipation 
Central  differences: 
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where  T  could  be  F*  or  Q  depending  on  the  form  of  the  equations 

Scalar  artificial  dissipation: 
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where  rj  is  even  and  one  more  than  the  order  of  accuracy 
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•  Runge-Kutta  Methods: 
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ESDIRK  Methods 


•  Explicit  first  stage  Singly-Diagonally 

Implicit  Runge-Kutta 

-  Stiffly  accurate 

-  Second-order  stage  accuracy 

-  FSAL  -  First  is  the  Same  As  Last 


ci  =  0 

0 

0 

0 

0 

0 

C2 

a  21 

A 

0 

0 

0 

C3 

&31 

• 

&32 

• 

A 

0 

0 

Ca-1 

• 

«(s-l)l 

• 

• 

a(s-l)2 

• 

• 

a(s-l)3 

A 

0 

cs  =  1 

h 

b2 

^3 

•  •  •  bs— l 

A 

h 

^2 

b3 

...  bs— i 

A 

A 

h 

A 

^2 

/\ 

^3 

A 

. . .  bg — i 

A 

bs 

Distribution  A  -  Approved  for  public  release;  Distribution  Unlimited 


7 


ESDIRK3  and  4 
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ESDIRK5 
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Von  Neumann  Analysis 

•  Often  used  to  study  stability  of  schemes 

•  Von  Neumann  analysis  is  used  to 
compare  schemes  for  accuracy 

-  Dissipation  error 

-  Dispersion  error 

•  Assumes  linear,  periodic  problems 

•  VNA  theory  and  more  results  are  in  the 
associated  paper 
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W  Dissipation,  CFL  =  1.0 
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Phase  Error 


Dispersion,  CFL  =  1.0 


Wave  Number 
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Magnitude  of  the  Growth  Factor,  |G| 


Dissipation,  CFL  =  10.0 


Wave  Number 
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Phase  Error 


Dispersion,  CFL  =  10.0 
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1  -D  Acoustic  Wave 


Unperturbed  Mach  number  of  0.5 

Poo  =  8.7077  x  lO-1^ 
puoo  =  1.7458  x  102  4r 
Too  =  400TT 
Too  =  2.871  X  102-^ 

7  =  1.4 

Perturbation  wave  -  20  points  per  wave  resolution 


Qo  Q OO  H-  M 5Qu,U±C 

A  A 

SQu,u±c  =  S'  COS  (foe) 

where  6  =  0.01 


•  More  results  in  the  paper 
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Amplitude  of  the  u+c  Characteristic  Variable 
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Amplitude  of  the  u+c  Characteristic  Variable 


1  -D,  CFL 


10.0, 1  Period 


Dissipation  Error 

Dispersion  Error 

Scheme 

VNA 

Simulation 

VNA 

Simulation 

Crank-Nicolson 

9.02  x  10“5 

2.44  x  10“a 

3.61  x  10 1 

3.61  x  10 1 

ESDIRK3 

4.99  x  10 1 

4.90  x  10 1 

1.92  x  10 1 

1.92  x  10 1 

ESDIRK4 

7.22  x  10-3 

7.25  x  10"3 

4.90  x  10"2 

4.90  x  10"2 

ESDIRK5 

5.10  x  10-2 

5.46  x  10-2 

1.38  x  10-2 

1.39  x  10-2 

Distribution  A  -  Approved  for  public  release;  Distribution  Unlimited 


17 


Amplitude  of  the  u+c  Characteristic  Variable 
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3-D  Isentropic  Vortex 


•  Free-stream  Mach  number  of  0.5 
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•  Perturbation  - 11  points  across  the  vortex 
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•  More  results  in  the  paper 
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Density 


3-D,  CFL  =  1.0,  40  Lengths, 
1 1  Points  Across  the  Vortex 
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Density 


3-D,  CFL  =  8.0 ,  40  Lengths, 
1 1  Points  Across  the  Vortex 
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Density 
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Conclusions 


•  2nd-  and  3rd- order  time  integrators  for  5th- order 
spatial  schemes  are  inadequate 

-  The  same  order  of  spatial  and  temporal  discretizations  is  preferable 

-  However,  ESDIRK5  is  not  much  better  than  ESDIRK4 

•  7  implicit  stages  vs.  5  implicit  stages 

•  Higher-order  time  integrators: 

-  Do  not  show  significant  improvement  on  coarse  grids  at  CFL  of  one 

-  Are  better  at  high  CFL  number 

-  Are  better  on  highly  refined  grids 

•  Spatial  error  usually  dominates  for  typical 
CFL  numbers  and  grid  resolutions 

-  Central  difference  plus  artificial  dissipation  schemes  are  inadequate 
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Future  Work 


•  Implement  more  accurate  spatial  schemes  of 
the  same  orders  of  accuracy 

-  Compact-difference  schemes 

-  Filtering  schemes 

•  Derive  better  ESDIRK  schemes  tailored  to  the 
desired  dissipation  and  dispersion  properties 

•  Add  preconditioning  to  take  maximum 
advantage  of  the  ESDIRK  time  integrators  for 
stiff  problems 

-  Improved  convergence  efficiency 

-  Improved  solution  accuracy 
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Questions??? 
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3-D,  CFL  =  8.0 
Different  Resolutions 
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